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SUMMARY 

Electrostatic (Langmuir) probes of both spherical and cylindrical geometry 
have been used to obtain electron number density and temperature In the exhaust 
of a laboratory arcjet. The arcjet thruster operated on nitrogen and hydrogen 
mixtures to simulate fully decomposed hydrazine In a vacuum environment with 
background pressures less than 0.05 Pa. 

The exhaust appears to be only slightly Ionized (less than 1 percent) with 
local plasma potentials near facility ground. The current-voltage characteris- 
tics of the probes Indicate a Maxwellian temperature distribution. Plume data 
are presented as a function of arcjet operating conditions and also position 
In the exhaust. 


INTRODUCTION 

An arcjet thruster Is an electrothermal propulsion device which heats a 
propellant by means of a high temperature arc discharge and then expands the 
hot gas through a nozzle to produce thrust. Because It offers relatively high 
specific Impulse at moderate thrust levels, the arcjet Is an attractive candi- 
date for auxiliary propulsion applications. At present, the power realistic- 
ally available for satellite propulsion requirements Is between 0.5 and 1.5 kW. 
When operated with storable propellants at these power levels, the arcjet can 
provide specific Impulses of 400 sec and higher, offering potentially large 
benefits for mission applications such as stationkeeping on communlcatldp sat- 
ellites. Recently, efforts have been Initiated, both at NASA and In Industry, 
to demonstrate the technology readiness of a low power (700 to 1500 W), dc arc- 
jet with storable propellants for such missions. Technology Issues currently 
being addressed Include arc start-up and stability, performance, lifetime, 
power processing, and spacecraft Integration (refs. 1 to 4). In conjunction 
with the above efforts, research programs to Investigate the plumes of arcjet 
thrusters have also begun. 

A detailed description of the arcjet exhaust Is required for not only a 
fundamental understanding of the plume character but also an accurate assess- 
ment of the potential Impacts that the plume may have on thruster performance 
and to the spacecraft upon Integration. The primary Integration concern Is the 
potential Impact due to EMI or Induced environment effects on communications. 
Thrust losses and thermal loading may also be of concern depending on the 
spacecraft/thruster configuration. The exhaust of an arcjet can be described 
as a low density, high temperature, supersonic flow field that Is slightly Ion- 
ized (less than a few percent) and chemically reactive In nature. Because of 
this complex combination of plume properties and also the fact that arcjet 
nozzles are typically very small. It Is difficult to obtain local, quantita- 
tive Information on the flow field with conventional probe techniques. 


This paper is declared a work of Ibe U.S. Government and is 
not subject to copyright protection in the United States. 


However, experimental Investigations of the arcjet exhaust have been con- 
ducted previously. In the early 1960's, Plasmadyne Corporation (ref. 5) meas- 
ured the shape and radiation Intensity of a 1 kW, hydrogen arcjet using motion 
picture photography and a thermopile radiometer. Additionally, McDonnell Air- 
craft Corporation (refs. 6 and 7) applied a variety of diagnostic techniques 
near the exit plane of high power (10 to 50 kW), hydrogen arcjets. The meas- 
urements Included: probe surveys to obtain electron number density and tem- 

perature, mass flux. Impact pressure, and stagnation enthalpy profiles; 
photometric techniques to measure exhaust velocity; and spectrophotometer 
observations to determine excitation temperature. Although a great deal of 
effort was expended In developing the above-mentioned techniques, much of the 
Information gained In the previous experiments was qualitative In nature and 
not directly transportable to present analyses. 

Lastly, Llmbaugh (ref. 8) determined the excited state density distribut- 
ion within the freely expanding plume of an argon arcjet (150 A, 20 V). He 

found the energy distribution did not fit a Maxwell -Boltzmann population of 

excited state number densities and that the electron temperature was much 
greater than the static gas temperature. Using spectral line absorption and 
emission measurements, electron number densities of 1 to 4xl0 14 /cm 3 were 
obtained at temperatures of approximately 0.6 eV. 

This paper presents Initial results In the experimental Investigation of 
low power, dc arcjet plumes. The primary research goals are: to obtain plasma 

number density and temperature as a function of position and operating condi- 
tion; to evaluate vacuum facility pressure effects on the measurements; and, 
ultimately, to assess the plume Impacts, subject to the above experimental 
results. For example, once the plume electron number density distribution Is 
known, estimates of the Impacts to communications signals transmitted through 
the plume can be calculated. The Langmuir probe Is the primary diagnostic 
technique to be used In this Investigation. 

The first portion of this paper provides a discussion of the theory and 
operation of Langmuir probes for measuring plasma properties. Details of probe 
design and fabrication, along with factors which govern the data analysis, will 
be provided. The experimental hardware and vacuum facility for the probe sur- 
veys are also discussed. 

The major portion of this paper presents Initial experimental data In the 
plume of a low power, dc arcjet. The arcjet used In this study Is a laboratory 
thruster, operating on a propellant mixture of nitrogen and hydrogen, to simu- 
late fully decomposed hydrazine. Data are presented as a function of position 
In the exhaust and also arc power. Comparisons are made between probes of dif- 
ferent size and geometry. Finally, recommendations for future work are 
proposed. 


NOMENCLATURE 


A probe surface area, m 2 

2 2 

A' wetted surface area of spherical probe (2*R p ), m 
A e nozzle exit area, m 2 
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A* nozzle throat area, m 2 

D probe diameter, m 

e electron charge, C 

G arc gap setting, m 

I current, A 

I 0 electron saturation current, A 

k Boltzmann constant, J/K 

L length, m 

1 cylindrical probe length, m 
m charged specie mass, kg 
M neutral mass, kg 

n number density, part1cles/m3 
R radius, m 

T temperature, K 

V voltage or potential, v 
x mean free path, m 

X-o Debye length, m 

<p nozzle half angle, degrees 
Sub- or Superscripts 
a arc 

c constrictor 

e electron 

f floating 

1 Ion 

n neutral particles 

p probe 

plasma 


00 


ANALYSIS 


The Langmuir probe Is one of the most widely used and accepted means of 
obtaining Information on plasma properties. In Its simplest form. It consists 
of a small conductive element and an associated power supply which Is capable 
of biasing the probe to various voltages, both positive and negative, with 
respect to the plasma to be Investigated. The current drawn to the probe as a 
function of voltage, or the "probe characteristic," yields two primary quanti- 
ties, electron number density and electron temperature. 

The theory that governs Interpretation of the probe characteristic Is 
quite complicated. Excellent reviews of general probe theory In both the pre- 
sence and absence of collisions may be found In references 9 to 11. 

Laframbolse (ref. 12) established both spherical and cylindrical probe theory 
In a colllslonless, stationary plasma. Further, the application of the 
Langmuir probe to flowing, colllslonless plasmas has been Investigated by 
French (ref. 13), Sonin (ref. 14), and Clayden (ref. 15). 

An Ideal probe characteristic Is Illustrated In figure 1. If the probe 
Is biased sufficiently negative with respect to the plasma potential, all elec- 
trons will be repelled and only Ions will be collected by the probe. A layer 
of positive charge, called an "Ion sheath," builds up around the probe. Its 
formation Is primarily dependent upon the Ion mass and temperature. As the 
probe voltage Is Increased positively, electrons are also attracted to the 
probe surface. At the floating potential, Vf, the magnitude of the Ion and 
electron currents are equal. Increasing the probe voltage further draws still 
more electrons to the probe. If the velocity distribution of electrons In the 
plasma Is Maxwellian, this transition region will be marked by an exponential 
Increase In electron current. At some point, V,,,, the probe potential Is equal 
to that of the plasma. There Is still a small Ion contribution at V«». Bias- 
ing the probe beyond V« causes electrons to be accelerated toward the probe 
and Ions to be repelled. In this region, the excess of negative charge near 
the probe surface again forms a sheath. As long as the probe size Is large In 
comparison with the Debye length, the sheath thickness may be assumed to be 
relatively constant with potential and equal to that of the probe dimensions. 


Provided that the probe does not disturb the plasma (or the macroscopic 
flow field) and that the electrons exhibit a Maxwellian distribution. Interpre- 
tation of the characteristic Is reasonably straightforward. In the transition 
region, the electron current collected by the probe will be given by 


A en 
P e 


/ kT 

Vss; 


exp 


- e < V n - V J 

P CO 

kT. 


( 1 ) 


Taking the logarithm of equation (1) and differentiating with respect to 
the probe voltage gives 


« < 1n v ■ it < 2 > 

e 

Consequently, the slope of a semi-log plot of electron current as a func- 
tion of probe voltage yields electron temperature. Plasma potential Is gener- 
ally obtained by extrapolating the linear portion of the transition and 
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electron accelerating regions of this same plot. At plasma potential, using a 
thin sheath analysis. 



(3) 


The saturation electron current, lot is the electron current at plasma potential 
(Vp » V,,,). Rearranging equation (3) gives the electron number density In terms 
of known quantities, as. 


n 


e 


(3.73X10 13 ) 


(4) 


The above analysis Is predicated on the assumption that 


X » Rp » \q 


(5) 


where the Debye length Is 



( 6 ) 


and x represents the collision mean free paths. A discussion of the relevant 
collision phenomena and the associated mean free paths may be found In 
appendix A. 


For low density plasma flows, previous Investigators (refs. 13 and 14) 
have used cylindrical probes aligned with the flow streamlines to avoid the 
effects of bulk fluid (neutral) motion. Spherical probes are promising because 
they do not have this alignment constraint. Clayden (ref. 15) has obtained 
reasonable results using spherical probes In a low density plasma jet; however, 
a "wetted* current collection area (one-half of the spherical surface area) was 
used In the calculation of n e (eq. (4)). This modification was made to 
account for the existence of a wake behind the probe which electrons are unable 
to penetrate. 


APPARATUS 
Langmuir Probes 

Langmuir probes of both spherical and cylindrical geometry were used In 
this Investigation. Figure 2 displays sketches of the stationary and movable 
probes. The stationary probes were made from carbide steel ball bearings with 
a radius of 0.562 cm. Each bearing was welded to a 0.76 mm tungsten wire. 
Stainless steel tubing was used to support the probe and also serve as an elec- 
trostatic shield. The ball bearing was located approximately 4 to 5 support 
diameters In front of the tube. Alumina tubing was used as the Insulator 
between the wire and the shield. It was recessed 0.16 cm from the end of the 
stainless shield. Kapton tape covered the shield In the vicinity of the 
sphere. A larger spherical probe (radius of 0.794 cm) was Installed at the 
off-axis stationary location for later measurements. 
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The support assembly for the movable probes had a U-shaped geometry 
because they had to sweep across the plume. The spherical probe was a carbide 
steel ball bearing of radius 0.316 cm. The cylindrical probe was made of tung- 
sten wire of radius 0.114 cm and length 0.952 cm. Both probes used ceramic 
Insulation as the main support at the location of the measurement. 

Figure 3 displays an electrical schematic of the probe bias circuitry. 

The bias supply could linearly ramp the probe voltage from a preselected nega- 
tive value to a positive value over a five second time period. The ramp volt- 
age ranges which could be selected were: -10 to +10 V, -50 to +50 V, and 

-100 to +100 V. The probe voltage was measured by a high Impedance divider 
network while the probe current was obtained by measuring the voltage across a 
1.00 Q resistor. The two signals were sent through Isolation amplifiers to 
an X-Y recorder. Because the current collected by the probe Increases exponen- 
tially near plasma potential, the probe bias supply ramped slowly enough so 
that the X-Y recorder could keep up. As long as the arcjet exhibited a stable 
operating condition, the dc current change with time varied by less than 
1 percent when the probe was biased at large positive or large negative volt- 
ages with respect to plasma potential. However, a slight difference In anode 
potential could have a marked effect In current collected near plasma poten- 
tial. Thermal and time dependent offsets In the amplifiers were zeroed out 
prior to each probe sweep. The estimated uncertainty In the current measure- 
ments was approximately 0.3 mA. 


Arcjet Thruster 

A photograph of the laboratory arcjet used In this Investigation Is shown 
In figure 4. The main body of the the thruster Is stainless steel with boron 
nitride Insulation. The thruster discharge region Incorporates a vortex sta- 
bilized, constricted arc design. The tungsten Insert shown In the figure 
serves both as anode and expansion nozzle. Major geometric characteristics of 
the arc constrictor and nozzle are found In table I. Reference 1 provides a 
complete description of the major design features of this thruster. 

Arcjet power was provided from a well regulated, pulse width modulated 
supply (ref. 2). The design consists of a starting circuit which provides 
brief, high voltage pulses to Initiate the arc discharge via Paschen breakdown 
and a second circuit for a well regulated, current-controlled mode of opera- 
tion. In normal operation, arcjet current ranged from 8 to 12 A, although most 
of the tests reported herein were conducted at 11 A. The arc current was meas- 
ured with a commercially available current probe while arc voltage was measured 
with a high Impedance divider network. The entire system was Isolated with 
anode potentials, during arcjet operation, typically near facility ground due 
to plasma bleed. The thruster floating voltage, anode to ground, was monitored 
with a 10 Hfi digital meter. 

Typically, the arcjet was operated on a 2:1 volume ratio of nitrogen to 
hydrogen to simulate fully decomposed hydrazine. The two gases were stored 
separately and mixed upstream of the thruster Inlet. Digital flow meters, cal- 
ibrated for the specific gas type, were used In conjunction with precision 
needle valves to control Individual volume flow rates. The full-scale reading 
of each digital meter was 5.00 slpm. The mixed propellant was fed to the 
thruster through a flexible feed tube. Unless otherwise noted, the mass flow 
rate of the thruster using the N 2 :H 2 mixture was 0.045 g/sec (2.00 SLPM N 2 
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and 4.00 SLPM H 2 ). Additionally, there were a number of Initial plume measure- 
ments made using only pure nitrogen as the propellant. The typical mass flow 
rate for these nitrogen-only tests was 0.061 g/sec. 


Instrumentation 

Measurements of the dc arc current and voltage, and the premixed propel- 
lant volume flow rates were monitored by a chart recorder. Isolation ampli- 
fiers and shielded cables were used for all signal processing. The real-time 
transients of the arc current were also monitored using an analog oscilloscope 
with a differential amplifier Input. Lastly, the thruster Inlet pressure, 
facility pressure, and thruster floating voltage were recorded manually at 
regular Intervals. Typical arcjet thruster operating conditions and perform- 
ance characteristics are listed In table II. 


Test Facility and Experimental Set-up 

The Langmuir probe surveys were conducted In a 1.5 m diameter by 5 m long 
vacuum facility at NASA Lewis Research Center. The tank Is equipped with four 
oil diffusion pumps and, for the the experimental conditions described herein, 
maintained a background pressure less than 5x10-2 p a (4x10-* Torr). 

The arcjet thruster, as displayed In figure 4, was mounted on a cradle 
assembly at the end of a push/pull rod that could be extended well Into the 
tank environment. Here, It rested on a fixed table to ensure that the thruster 
position would be stable and repeatable from one test date to the next. The 
thruster exhausted along the tank centerline as shown In figure 5. The two 
stationary, spherical probes were mounted Inside the tank at the locations 
denoted In the figure. One probe was located on the plume centerline, 30.5 cm 
downstream of the thruster exit plane, while the other probe was at the same 
axial plane but 15.2 cm off centerline. The movable probes were positioned In 
the plume using a push/pull rod entering through a side port entrance. Probes 
of both spherical and cylindrical geometry were used. The only rotation avail- 
able with this arrangement was In the vertical plane; therefore, cylindrical 
probe measurements could only be obtained on the plume centerline. The spheri- 
cal probes were used to obtain radial profiles of the arcjet plume properties. 
The axial distance for these measurements was 18.4 cm. 


PROCEDURE 

The positions of the thruster and the probes were fixed at their predeter- 
mined locations prior to the experiment and every time a probe was changed. 

The Langmuir probes were cleaned chemically before placement In the tank. Sub- 
sequent cleaning at the beginning of each test period was accomplished by Ion 
sputtering. The probes were biased to a large negative potential for several 
minutes while the thruster was operating. Generally, the thruster was allowed 
to run for 20 to 30 min after the arc discharge was Initiated so that It could 
reach a steady-state condition. Subsequent changes to thruster operation, such 
as variations In arc current, required only a few minutes for the arcjet to 
reach a new steady-state operating condition. 
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Nitrogen-Hydrogen Mixture Operation 

The bulk of data presented were taken with the arcjet operating on a 2:1 
mixture ratio of hydrogen and nitrogen, as previously specified. The Langmuir 
probe drive supply was tied to anode potential and the probes were ramped line- 
arly through a voltage range of -10 to +10 V. Since the greatest change In 
current collected by the probe varies over a short voltage range, probe data 
were taken with the X-Y recorder at two different sensitivities. 

Using the stationary probes, probe data were obtained at the two down- 
stream locations as a function of arcjet operating condition. With the mass 
flow constant at 0.045 g/sec, the arc current was varied from 9 to 12 A. Typi- 
cal arc voltages for these measurements ranged from 100 to 105 V, although 
sometimes excursions between 92 and 110 V were noted for a given arc current 
and a fixed flow rate. Since the arc voltage Is primarily dependent upon the 
manner In which arc attachment Is made at the anode, slight variations In arc 
voltage are not unusual. As long as the arcjet exhibited a stable operating 
condition, with a stable arc voltage and a nonwavering plume, plume data could 
still be obtained. Slight voltage excursions do not significantly Impact 
thruster performance (ref. 16). 

Plume data were often taken under similar arcjet operating conditions on 
different test dates and at different times during a given test date to verify 
repeatability and consistency of the measurements. Also, data were taken with 
probes of different size and geometry. In all cases, the arcjet operating con- 
dition was established by setting the arc current (usually 11 A) and mixture 
flow rate (0.045 g/sec) to the same values. 

A radial profile was also obtained with a movable probe of spherical geom- 
etry. It took a period of approximately 15 min to obtain all of the radial 
data. During that time the arcjet remained stable In terms of Its operating 
characteristics, Including arc voltage. Nonetheless, stationary probe data 
were obtained before and after the the radial test sequence to ensure consis- 
tent plume characteristics. 


Nitrogen Only Tests 

A limited amount of experimental data was obtained with the arcjet operat- 
ing on pure nitrogen gas. The arc current for the data presented was 10 A at 
a flow rate of 0.061 g/sec. The Langmuir probes were tied to cathode common, 
which typically ran at a potential of 60 V below facility ground. The probes 
ramped through a -100 to +100 V sweep at two different recorder sensitivities. 


Experimental Data Reduction 

A computer code originally developed at the Colorado State University 
(ref. 17) was used to calculate the plume parameters from the V-I characteris- 
tic. Probe electron current (total probe current minus Ion contribution) and 
associated voltage pairs were Input In the program. A least squares fit was 
used to extrapolate the linear portions of the retarding and accelerating 
regions of the probe characteristic. Plasma potential was defined as the point 
at which the two extrapolated lines Intersected. The computer code relied on 
an estimate of plasma potential and also the voltage at which the accelerating 
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portion of the characteristic began. Consequently, the data were first plotted 
graphically to obtain reasonable estimates of these values. The plasma poten- 
tial Input was chosen as the value of probe voltage at which the slope of the 
logarithmic current versus voltage first deviated from linearity. 


RESULTS AND DISCUSSION 

The first portion of this section discusses typical Langmuir probe charac- 
teristics for arcjet operation on pure nitrogen and also nitrogen-hydrogen mix- 
tures. The general character of the plume under the two different modes of 
operation Is described, along with estimates of the uncertainty of the measure- 
ments. Comparisons are made between probes of different size and geometry. 
Lastly, the exhaust of an arcjet operating on mixtures Is described In terras 
of electron number density and temperature as a function of arcjet operating 
condition and also location In the plume. 


Representative Probe Characteristics 

Figures 6 and 7 display representative V-I characteristics for the spheri- 
cal probes located 30.5 cm downstream of the thruster exit on the plume center- 
line. In the two figures, the logarithm of the current collected by the probe 
(total probe current minus Ion contribution) Is plotted against the probe volt- 
age with respect to arcjet common. The different symbols In the figures denote 
the X-Y recorder sensitivities. 

In figure 6, the arcjet was operating on pure nitrogen propellant at a 
mass flow rate of 0.06 g/sec. Thenarc current was 10 A with a corresponding 
arc voltage of 63 V. The anode potential Is approximately 3 V above facility 
ground. 

The linearity over three decades of current change Indicates a Maxwellian 
distribution of electron velocities. Applying the foregoing analysis, extrapo- 
lation of the linear portions of the transition and accelerating regions' yields 
an electron number density of 1.9xl0 10 /cm 3 at a temperature of 0.86 eV. 

The estimated uncertainty In the number density calculation Is about 8 percent, 
primarily due to the uncertainty In the slope determination (±0.10 eV). The 
plasma potential appears relatively close to that of anode potential, at 2 V. 

Similarly, figure 7 displays the characteristic for the same probe with 
the arcjet operating on nitrogen-hydrogen mixtures. In this figure, the probe 
voltage Is referenced to the arcjet anode, which Is again fairly close to 
facility ground, at 4-0.8 V. 

Even though the arc power Is considerably higher In this case, the plume 
exhibits a much lower concentration of charged species. Because the current 
collected by the probe Is linear over at least two decades, a Maxwellian dis- 
tribution Is assumed, and the electron density Is calculated at 3.5xlOVcm 3 . 
However, the estimated uncertainty Is 14 percent, since overall current values 
are considerably smaller than In the previous case. Also, the plume electron 
temperature Is slightly cooler with the arcjet operating on mixtures. The 
measured electron temperature Is 0.51 eV ± 0.10 eV. The Oebye length, based on 
the measured values of T e and n e . Is 9.0xl0~ 3 cm, a factor of sixty smaller 
than the probe radius, so the thin sheath analysis Is applicable. Values of 
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the charged-charged and charged-neutral mean free paths are at least one to two 
orders of magnitude larger than the probe radius at this measurement location. 
Estimates of the relevant mean free paths for the various collision phenomena 
may be found In table A1 . 

Variation In electron number density In the plume Is due to the combined 
effect of gas dynamic expansion and also recombination. Recombination may be 
the result of radiative, dissociative, or three body collisions In the exhaust. 
The reason for the lower concentration of electrons In the plume of the arcjet 
operating on mixtures Is presumably due to the large amount of hydrogen In the 
flow. Hydrogen, although It has the same Ionization potential as nitrogen, has 
significantly higher electron-ion recombination coefficients. Consequently, 
although the collision frequencies In the plume of the arcjet operating of pure 
nitrogen are probably higher than In the case of mixtures, as there are higher 
neutral densities, a smaller percentage of the charged species recombine. The 
data obtained thus far Is Insufficient to fully document the relative magni- 
tudes of colllslonal effects on recombination In the arcjet exhaust. 


Plume Data With Arcjet Operation on Mixtures 

All of the experimental data presented herein apply to arcjet operation 
with nitrogen and hydrogen mixtures to simulate decomposed hydrazine. Figure 8 
displays electron number density values obtained with the stationary spherical 
probes. The data In figure 8(a) were taken on the plume centerline, 30.5 cm 

downstream of the thruster exit, while the data In figure 8(b) were taken at 

the same axial station but 15.2 cm off axis. Because a set arc current did not 
always give a consistent arc voltage, the calculated electron number density 
values are plotted as a function of arc power. The various symbols, as 
described In the figures, denote the different arc current settings. The esti- 
mated uncertainty In the data Is 14 percent. Both figures display an Increase 

In measured electron number density with arc power. Further, within the limi- 
tations of the probe technique, there Is a reasonable consistency In the meas- 
ured values of electron number density for a given arc current and voltage. 

For example, the closed symbol data taken at 10 A (fig. 8(a)) were taken during 

the same test, but at time Intervals of at least 40 to 80 min apart. Further, 
the open symbol data point at 10 A was taken during a different test date but 
at approximately the same operating condition. 

Figure 9 displays corresponding values of electron temperature for the 
stationary probe on the plume centerline. A systematic variation In tempera- 
ture with arc current, voltage, or power Is not evident, as the scatter In the 
data Is fairly large. One might expect an Increase In electron temperature 
with arc current, because the core temperature Increases with current but, 
within the scatter of data found, this trend was not directly apparent, as 
shown figure 9. However, It should be noted that the average temperature for 
measurements at 10 A Is 0.53 + 0.10 eV, while at 11 A, the average value Is 

slightly higher at 0.64 + 0.10 eV. 


Variations In Probe Size and Geometry 

Plume data were taken with spherical probes of two different radii at the 
off axis, stationary location. Referring back to figure 8(b), there Is excel- 
lent agreement In measured number density between the two probes under similar 
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arcjet operating conditions. However, differences do exist In the measured 
electron temperatures. Table III provides specific Information on Individual 
data points. 

A comparison of probe geometry was also made using the movable probe 
arrangement. Table IV presents Individual probe data taken with the spherical 
and cylindrical probes at an axial distance of 18.4 cm, on the plume center- 
line. The electron temperatures are relatively consistent Independent of the 
probe geometry but there Is a definite Increase In measured electron nunb«r 
densities for the spherical probe data. There Is an Increased uncertainty In 
the cylindrical probe data due to the very small amount of current drawn to the 
probe with bias voltage, but this In Itself does not account for the discrep- 
ancy. Clayden (ref. 15) also found that spherical probes gave higher values 
of electron number densities than cylindrical shapes. It should be mentioned 
that a wetted surface area (one half of the spherical surface area) was used 
to Interpret the spherical probe data. While this modification Is made to 
account for the existence of a wake behind the probe which electrons are unable 
to penetrate (ref. 15), It may represent an a source of error In Interpretation 
of data and certainly warrants further Investigation. 


Radial Plume Measurements 

A spherical probe was used to obtain a radial profile at an axial distance 
of 18.4 cm downstream of the thruster exit plane. Figure 10 presents the meas- 
ured electron number density and temperature variation as a function of angle 
off plume centerline. As shown In figure 10(a), the measured number density 
decreases exponentially. The measured number density Is approximately one half 
of the centerline value at an angle of 25°. 

Figure 10(b) presents electron temperature values based on the same set 
of measurements. There Is no apparent variation In temperature within the 
radial distance off centerline (~23 cm maximum). These findings are similar 
to previous Investigations (refs. 7 and 15). Even though the neutral gas tem- 
perature varies markedly with radial position, the electrons do not follow this 
pattern because of the poor energy exchange between the heavy particle neutrals 
and the highly mobile electrons (ref. 14). 


CONCLUDING REMARKS 

Langmuir probes of spherical and cylindrical geometry have been used to 
Investigate the exhaust of a low power, dc arcjet thruster operating on nitro- 
gen and hydrogen mixtures. Preliminary measurements with a spherical probe 
located ~30 cm downstream of the thruster exit plane have Indicated electron 
number densities on the order of 3xl0 9 /cm 3 at a temperature of 0.5 eV. 

This low concentration of electrons represents less than 1 percent of the total 
centerline density. An Increase In the measured number density was observed 
with Increased arc power. 

A radial profile at a distance of 18 cm downstream of the thruster was 
also obtained. The measured electron density profile decreased exponentially 
as a function of angle off centerline while the values of electron temperature 
exhibited no systematic variation with radial position. Probes of different 


11 


size yielded consistent results, within the limitations of the existing meas- 
urement technique. However, a spherical probe gave higher measured electron 
densities than did a cylindrical shape, although the measured electron tempera- 
tures agreed. 

It Is apparent that detailed plume surveys of the arcjet exhaust will 
require probes of different size and geometry. In some regions It may be 
Impossible to avoid the use of complex probe theories In order to Interpret the 
probe characteristics. The macroscopic effects of the flow field on sheath 
formation (thick sheath analysis) and also the application of continuum methods 
require Investigation. Further, although the spherical probe does not require 
alignment with flow streamlines, there Is some question about the true current 
collection surface area near plasma potential. 

To obtain better measurement of Ion saturation current, It will be neces- 
sary to obtain greater accuracy In the current collected by the probe, probably 
down to the tens of micro-amperes. This will permit Information on Ion temper- 
atures and densities to also be extracted from a single probe characteristic. 
Lastly, vacuum facility pressure effects play a large role In the region where 
meaningful Information can be obtained. The effects of facility pressure on 
the gas dynamic expansion and also the recombination rates will require 
detailed Investigation. 
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APPENDIX A - COLLISION PHENOMENA 

In order to apply collisionless theory to the Langmuir probe characteris- 
tic, It Is Important that both the probe radius and the Debye length are small 
In comparison with the representative mean free paths. In a plasma where the 
percentage of Ionization Is small, l.e., less than a few percent, collisions 
between charged species and also neutrals become relevant. Tables A1 and A2 
list order of magnitude estimates of the most Important mean free paths for the 
two primary plume centerline measurement locations. The collisions between 
charged species were calculated using Spltzer's relation (refs. 14 and 18) 
while the mean free paths between neutrals were calculated using gas kinetic 
theory. Electron-neutral collisions (both elastic and Inelastic) were calcu- 
lated using Maxwellian averaged cross sections (ref. 19) based on emperlcal 
data for the collision cross sections as a function of electron velocity. Ion- 
neutral mean free paths were estimated using the assumption that x^ n Is 
approximately one-third of \ nn (ref. 14). 

The available experimental data were used to obtain necessary estimates 
of charged specie number densities at the two measurement locations. Reason- 
able values of the neutral densities, of both molecular and dissociated 
species, were not easily attainable. A detailed plume analysis requires a com- 
plete description of the nozzle flow field as Input. The highly viscous arcjet 
nozzle flow, coupled with the unwieldy size of the nozzles, precludes the 
Immediate adaptation of conventional nozzle flow computational techniques. 
Further, the nonequilibrium nature of the flow, In conjunction with the fact 
that a considerable amount of heat addition occurs In the nozzle Itself, add 
considerable complexity to the problem. 

The following technique was employed to describe the expansion of the neu- 
tral gas mixture Into vacuum. A source flow analysis (ref. 20) which has dem- 
onstrated reasonable agreement with measurements of similar flow fields was 
adapted. The technique Is a modified version of Simons' method which describes 
the density at a given location In the plume In terms of a centerline density 
which Is calculated from continuity relations. Major plume parameters are 
written In terms of conditions at the nozzle exit and also the boundary layer 
thickness. 

A 1-D, Isentroplc analysis of the nozzle flow would be Invalid In this 
situation. The arcjet nozzle temperature profile Is very steep along the noz- 
zle centerline; however, the majority of the gas flow passes through the rela- 
tively cooler boundary layer. The centerline Mach number will probably be 
reduced and the viscous nature of the flow will probably lend Itself to a very 
large portion of the boundary layer as being subsonic. For this analysis, 
experimental performance data were used to obtain nozzle exit properties. The 
boundary layer thickness was assumed to be ~70 percent of the exit radius. 
Conditions at the nozzle throat were calculated using a discharge coefficient 
of 0.9 and an assumed Mach number of 1. 
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TABLE I. - ARC JE1 THRUSTER DIMENSIONS 


Nozzle geometry 

Exit radius, m R e = 4.76x10-3 

Throat radius, m R* = 3.18x10-* 

Half angle, deg tp = 20 

Area ratio A e / A * = 225 

Arc constrictor geometry 

Constrictor diameter, m D c * 6.35x10-* 

Constrictor length, m L c = 6.35x10-* 

Arc gap setting, m 6 = 7.37x10-* 


TABLE II. - TYPICAL ARCJE1 THRUSTER PERFORMANCE CHARACTERISTICS 


Hp:N 2 mixture ratio 

2:1 

2:1 

3:1 

3:1 

N2 flow rate, slpm 

2.00 

2.00 

1.70 

1.70 

H 2 flow rate, slpm 

4.00 

4.00 

5.10 

5.10 

Total mass flow rate, g/sec 

0.044 

0.044 

0.041 

0.041 

Average molecular weight 

10.6 

10.6 

8.5 

8.5 

Inlet pressure, N/m? 

4.58x1 05 

4.66x105 

4.53x105 

4.71x105 

Arc voltage, V 

93.0 

95.4 

95.3 

92.5 

Arc current, A 

10.0 

11 

10 

12 

Arc power, W 

930 

1049 

953 

1110 

Thrust, g 

17.3 

18.0 

17.0 

17.6 

Specific Impulse, sec 

393 

409 

414 

430 

Overall efficiency, percent 

34 

33 

34 

32 

Cold flow Isp, sec 

114 

114 

125 

125 

Cold flow thrust, g 

5.0 

5.0 

5.1 

5.1 
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TABLE III. - COMPARISON OF LANGMUIR PROBE SIZE VARIATION 


[Location: x = 30.5 cm; R = 15.2 cm.] 


Thruster operating condition 





Arc current, A 

10.9 

10.9 

10.8 

10.9 

Arc voltage, V 

100 

104 

104 

103 

Arc power, W 

1090 

1134 

1123 

1123 

Mass flow rate, slpm 

2. 0/4.0 

2. 0/4.0 

2. 0/4.0 

2. 0/4.0 

Mass flow rate, g/sec 

0.045 

0.045 

0.045 

0.045 

Inlet pressure, N/m 2 

4.8xT0 5 

4.8xl0 5 

4.8x105 

4.8xl0 5 

Tank pressure, N/m 2 

4.9x10*2 

5.1x10-2 

5.2x10-2 

5.2x10-2 

Floating voltage, V 

*0.4 

*0.1 

*0.3 

+0.4 

[(<•) to grd] 

Langmuir probe Information 





Probe type 

Sph 

Sph 

Sph 

Sph 

Probe radius, cm 

0.562 

0.562 

0.794 

0.794 

Plume characteristics 





Electron temperature, K 

1 .2x10^ 

1 .2x104 

6.3x103 

7 . 7x1 0 3 

Electron temperature, eV 

1.0 

1.1 

0.54 

0.66 

Electron number density, cm- 3 

1 .7x109 

2.0x109 

2.0x109 

1 .6x109 

Plasma potential, V [with 

+ 1 .4 

+ 1 .1 

+0.9 

*1 .0 

respect to facility grd] 





Electron sat. current, mA 

9.0 

10.7 

15.6 

13.4 

Debye length, cm 

1 .8xl0- 2 

1.7x10-2 

1 .2x10-2 

1 .5x10-2 

Rp/\ D 

31 

33 

66 

53 

Xln/Rp 






TABLE IV. - COMPARISON OF LANGMUIR PROBE GEOMETRY VARIATION 


[Location: x = 18.4 cm; plume centerline.] 


Thruster operating condition 







Arc current, A 

10.1 

10.9 

10.9 

10.9 

10.9 

10.8 

Arc voltage, V 

108 

98 

103 

94 

94 

92 

Arc power, W 

1091 

1068 

1123 

1025 

1025 

994 

Mass flow rate, slpm 

2. 0/4.0 

2. 0/4.0 

2. 0/4.0 

2. 0/4.0 

2. 0/4.0 

2. 0/4.0 

Mass flow rate, g/sec 

0.045 

0.045 

0.045 

0.045 

0.045 

0.045 

Inlet pressure, N/m 2 

4.7x105 

4.8x105 

4.9x105 

4.8x105 

4.8x105 

4.7x105 

Tank pressure, N/m 2 

5.0x10*2 

4.9x10-2 

4.9x10-2 

5.3x10-2 

5.2x102 

5.1X10 2 

Floating voltage, V 

4-0.6 

*0.4 

+0.4 

*0.1 

0.0 

0.0 

[(+) to grd] 

Langmuir probe Information 







Probe type 

Sph 

Sph 

Sph 

cyl 

cyl 

cyl 

Probe radius, cm 

0.316 

0.316 

0.316 

0.057 

0.057 

0.057 

Probe length, cm 




0.476 

0.476 

0.476 

Plume characteristics 






1 .5x104 

Electron temperature, K 

8.6X10 3 

1 .0x1 0* 

8.0x103 

1 .2x104 

8.5x103 

Electron temperature, eV 

0.74 

0.86 

0.69 

1.0 

0.73 

1 .3 

Electron number density, cm- 3 

8.0x109 

6.2x109 

7.8xl0 9 

2.5x109 

3.0x109 

2.6x109 

Plasma potential, V [with 

*1.5 

+0.4 

+1 .4 

*1 .0 

*0.8 

*0.7 

respect to facility grd] 
Electron sat. current, mA 

11.6 

9.5 

10.8 

2.3 

2.4 

1 .4 

Debye length, cm 

7.1xl0- 3 

8.8x10-3 

7.0x10-3 

1 .5xl0- 2 

1 .2x10-2 

1 .7x10-2 

Rp/Xp 

45 

36 

45 

4 

6 

3 

Xln/ R p 

63 ; 

63 

63 

350 

350 

350 
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TABLE Al. - RELEVANT COLLISION PHENOMENA AND ESTIMATED MEAN FREE PATHS 
[Location: Plume centerline, X » 18.4 cm.] 


Type of 

Collision 

Mean free 

Assumptions 

collision 

partners 

path, 

cm 



Elastic 


Ion-Ion 

electron-electron 

electron-ion 

Ion-electron 


electron-N 2 5xl0 2 

electron-N 5x1 0 3 

electron-H 2 IxlO 2 

electron-H 3xl0 2 

1on-N 2 3X10 1 

Ion- N 2xl0 2 

1on-H 2 2x1 0 1 

Ion- H 1 xlO 2 

N 2 -N 2 9x1 0^ 

N 2 -H 2 6x1 0^ 

N-N 7x101 

H-H 5x102 

Inelastic electron-N 2 (lonlz.) 3x1 05 

electron-H 2 

electron-N 2 (exclt.) 3x10® 
electron-Ho 


use rms v^; kT^ = 0.1 eV; n^ = 8xl0 9 
use rms ve; kT e = 0.8 eV; n e = 8xl0 9 

*»e1 ~*-ee ( re ^ • 14) 
n 1 = 8x1 0 9 (refs. 18 and 14) 
kT e * 0.8 eV 

Q * 7x10“® cm 3 /sec; n N = 2.1xlQl2 cnr 3 
2 

Q = 9x1 0~® cm 3 /sec; n N * 4.6xlOH cnr 3 

Q = 2x10“® cm 3 /sec; n H = 4.0xlQl2 cm- 3 
2 

Q * lxl0~ 7 cm 3 /sec; n H = 1 . 4x1 0 3 2 cm' 3 
Mn ^0.3 X nn 


nij = 2.1xl0l2/cm 3 
2 

nu = 4.6x101 2/cm 3 
2 

n N = 4.6xlOH/cm 3 
n^ = 1. 2x1 0l 2/cm 3 
Q = 7x10-13 cm 3 /sec 


Q = 4x10-12 cm 3 /sec 


TABLE A2. - RELEVANT COLLISION PHENOMENA ANO ESTIMATEO MEAN FREE PATHS 
[Location: Plume centerline, X = 30.5 cm.] 


Type of 
collision 

Collision 

partners 

Mean free 
path, 
cm 

Assumptions 

Elastic 

Ion-Ion 

4X10 1 

use rms v^; kT^ = 0.1 eV; = 4xl0 9 /cm 3 


electron-electron 

1x1 0 2 

use rms ve; kT e = 0.8 eV; n e - 4xl0 9 /cm 3 


electron- Ion 

IxlO 2 

*-e1 ~*ee ( r ©f- 14) 


Ion-electron 

8xl0 3 

use rms v^; n^ = 4xl0 9 /cm 3 




kT e * 0.5 eV 


electron-N 2 

IxlO 3 

0 = 5x1 O'* 8 cm 3 /sec; = 9x1 0 11 /cm 3 


electron- N 

2x10* 

Q = IxlO' 8 cm 3 /sec; n^ = 2xl0H/cm 3 


electron~H 2 

4xl0 2 

Q = 8x!0" 8 cm 3 /sec; nn * 1 .7xl0 12 /ctn 3 


electron-H 

7xl0 2 

Q = IxlO' 7 cm 3 /sec; n^ = 6xlOH/cm 3 


1on-N 2 

7x10^ 




lon-N 

4xl0 2 






r Mn '‘■'0-3 *-nn 


1on-H 2 

6x10^ 




lon-H 

2xl0 2 > 




N 2 -N 2 

2xl0 2 

nw = 9xl0H/cm 3 
2 


n 2 -h 2 

4xl0l 

n^ - 9xl0ll/cm 3 ; n H ^ = 1 . 7x1 Ol 2 / C m 3 


h 2 ~h 2 

2xl0 2 

n^ = 1 .7x10^ 2 /cm3 


N-N 

1 xlO 3 

2 

n# = 2x1 0 11 /cm 3 


H-H 

6x1 0 2 

* 6x1 Ol 1 /cm 3 

Inelastic 

electron-N 2 ( lonlz. ) 

5x1 05 

Q = 6x1 O'" 1 3 cm 3 /sec 


electron-H 2 




electron-N 2 (exclt.) 

7xl0& 

Q = 4xl0~l 4 cm 3 /sec 


electron- H 2 








FIGURE 1. - IDEAL LANGMUIR PROBE CHARACTERISTIC WITH RESPECT TO AN 
ARBITRARY REFERENCE POINT. 
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BALL BEARING, | 
0,562 CM RADIUS- 1 


(A) STATIONARY LANGMUIR PROBES. 



TUNGSTEN WIRE. 
0.057 CM RADIUS 

0.952 CM LENGTH 7 

/ 



(B) MOVABLE LANGMUIR PROBES OF SPHERICAL AND CYLINDRICAL GEOMETRY. 
FIGURE 2. - ILLUSTRATIONS OF LANGMUIR PROBE TYPES. 
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FIGURE 3. - ELECTRICAL SCHEMATIC OF LANGMOIR PROBE DRIVE CIRCUIT 
AND DATA ACQUISITION. 
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FIGURE 1. - PHOTOGRAPH OF ARC JET THRUSTER MOUNTED ON CRADLE ASSEMBLY. 
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FIGURE 5. - SKETCH OF EXPERIMENTAL SET-UP (TOP VIEW). 
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PROBE VOLTAGE WRT <-), V 

FIGURE 6. - TYPICAL V-I CHARACTERISTIC FOR SPHERICAL PROBE; 
x = 30.5 CM/ PLUME CENTERLINE; ARCJET-N 2 PROPELLANT. 
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FIGURE 7. - TYPICAL V-I CHARACTERISTIC FOR SPHERICAL 
PROBE; X - 30.5 CM, PLUME CENTERLINE; MIXTURE. 
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ARC POWER, W 

(B) LOCATION: x = 30.5 CM, R = 15.2 CM. 


FIGURE 8. - ELECTRON NUMBER DENSITY AS A FUNCTION OF ARC 
POWER. PROPELLANT; N 2 -H 2 MIXTURE. 
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ARC POWER, W 


FIGURE 9. - ELECTRON TEMPERATURE AS A FUNCTION OF ARC POWER. 
LOCATION: X = 30.5 CM, PLUME CENTERLINE. PROPELLANT: 
N 2 -H 2 MIXTURE . PROBE RADIUS: 0.562 CM. 
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ANGLE OFF CENTERLINE, deg 

(B) ELECTRON TEMPERATURE AS A FUNCTION OF ANGLE OFF 
CENTERLINE. 

FIGURE 10. - RADIAL PLUME SURVEY AT X = 18.4 CM. PROPEL- 
LANT: N 2 -H 2 MIXTURE; PROBE RADIUS = 0.316 CM. 
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